Non-oscillating solutions to uncoupled Ermakov 
systems and the semiclassical limit 



A. Matzkin 

Abstract. The amplitude-phase formulation of the Schrodinger equation is 
investigated within the context of uncoupled Ermakov systems, whereby the amplitude 
function is given by the auxiliary nonlinear equation. The classical limit of the 
amplitude and phase functions is analyzed by setting up a semiclassical Ermakov 
system. In this limit, it is shown that classical quantities, such as the classical 
probability amplitude and the reduced action, are obtained only when the semiclassical 
amplitude and the accumulated phase are non-oscillating functions respectively of the 
space and energy variables. Conversely, among the infinitely many arbitrary exact 
quantum amplitude and phase functions corresponding to a given wavefunction, only 
the non-oscillating ones yield classical quantities in the limit ?i — > 0. 
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1. Introduction 

Systems of the form 

dMt) + k\t)n{t) = -^Y{a{t)/u{t)) (1) 

dMt) + k\t)a{t) = -^Z{uit)/ait)), (2) 

where Y and Z are arbitrary functions of their arguments, are generically known as 
Ermakov systems. They are characterized by the existence of a first integral, the 
Ermakov (or Lewis- Ray- Reid) invariant hnking the solutions of Eqs. (|l|) and (|]), thereby 
giving rise to the so-called nonlinear superposition principle. Intensive studies of their 
properties, such as their linearization or their generalization to higher dimensions 
0, ^ have been undertaken to extend the remarkable results concerning uncoupled 
systems, ie Y{Q = and Z^Q = a^C, where a is a constant, obtained by Ray and Reid 

Although the paradigm in physical applications of uncoupled Ermakov systems has 
been the classical linear time-dependent harmonic oscillator (t being the time variable), 
where the full power of Hamiltonian structure 0, Lagrangian mechanics and Noether 
symmetries have been employed, it has also been remarked that uncoupled systems 
link the time-independent linear Schrodinger equation to a nonlinear 'auxiliary' equation 
in the following way: 

h^dlu{x) +p'^{x)u{x) ={) (3) 
fi^d'la{x) +p^{x)a{x) = ^^^^^ , (4) 

where x refers to the space variable, and p{x, E) to its conjugate momentum in classical 
mechanics. E is the energy, assumed to be conserved, of the system. We have recently 
shown that the non-linear equation (^) corresponds to the equation for the amphtude 
function a{x) in the amplitude-phase formulation of the Schrodinger equation, which 
arises by performing a so-called Milne transform on the wave-function; the phase 
function (f){x) is obtained by integrating the relation dx4> = (see Sec. 2). 

We shall be concerned throughout this paper by the oscillatory properties of the 
solutions of the uncoupled Ermakov system formed by Eqs. (^ and (^ in the specific 
case of a potential energy function having a single minimum. Actually, since Eq. (|^) 
defines a Sturm-Liouville problem, the oscillatory properties of the solutions are well- 
known , and we will restrict our analysis to the oscillations as a function of x and E, 
of the amplitude and phase functions. Our first aim will be to show that, absolutely 
smooth, that is non-oscillating amplitude-phase functions can be constructed. Let us 
recall that the amplitude-phase formulation of the one-dimensional Schrodinger equation 
is frequently used in quantum scattering theories that explicitly include closed channels, 
such as quantum defect theory (see and references therein). In those situations, it is of 
prime importance for a{x, E) to be a smooth function of both the space (usually radial) 
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coordinate x and the energy E, since the scattering parameters (for example the phase- 
shifts) are defined in terms of amphtude-phase functions. However, by the principle 
of nonlinear superposition (Sec. 2.1), a may be expressed in terms of 2 independent 
solutions Ml and U2 of Eq. (^) - solutions which as known oscillate between the turning 
points of the potential. It follows that a{x) generally oscillates between the turning 
points. In practical implementations of amplitude-phase formalisms, numerical methods 
aiming at minimizing the amplitude of the oscillations have been devised. We have 
proposed in [§] such a method based on the invariant and the nonlinear superposition 
principle in the context of Ermakov systems. Here, the main point to be examined 
consists in the relationship between non-oscillating amplitude-phase functions and the 
functions obtained in the semiclassical {h —>■ 0) limit. More specifically, we will prove 
that in this limit, the only non-oscillating solutions are the ones that yield classical 
quantities: in particular, the only semiclassical phase function that does not oscillate is 
the classical reduced action, and conversely the quantum continuation, for finite h, of 
the classical reduced action is a non-oscillating function. 

To this end, some properties of amphtude-phase functions, their behaviour as a 
function of x and E, as well as the connection with Ermakov systems will be recalled in 
Sec. 2. In Sec. 3, we shall discuss semiclassical amplitude-phase functions by setting up 
a semiclassical Ermakov system; the classical probability amplitude and reduced action 
will appear as a particular solution of the semiclassical amplitude and phase functions. 
Those results will then be employed in a formal h expansion of the quantum amplitude- 
phase functions, to prove that provided the first order functions are non-oscillating, the 
solutions to each order in h will then not oscillate (Sec. 4). This will be followed in Sec. 
5 by a discussion of the results, in particular in relation to quantization of classically 
integrable systems. 

2. Amplitude-phase functions 

Our concern here will be the 'auxiliary' amplitude and phase functions of the Schrodinger 
Eq. d^), with p^(x) = 2m{E — V{x)), ie a particle of mass m and energy E trapped in 
a potential well V{x) having a single minimum. V{x) is defined on an interval ]si,S2[ 
(typically Si^2 = ±C)0 or 0). Atomic units will be used throughout, except for the h 
factors which will be reestablished where appropriate. 

2.1. Nonlinear superposition principle 

We collect in this paragraph the main results concerning amplitude-phase functions 
that will be useful in what follows, omitting details (for more details and the relevant 
references, see [§]). Eq. (^ defines a Sturm-Liouville problem, typically a singular 
problem on the half line or real line when vanishing boundary conditions at Si and 
S2 are implemented. However, our interest lies not in the specific eigenf unctions or 
eigenvalues of the Schrodinger equation, but in relating linearly independent solutions 
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of the hnear equation to amphtude and phase functions. We will denote by Ui and 
U2 two independent solutions of Eq. (^) respectively regular at si and S2 and with 
Wronskian W = >V['Ui,-U2] = {dxUi)u2 — ui{dxU2) ■ 

A general solution u{x) of Eq. which is readily written in terms of the 

independent solutions ui and U2, may also be obtained as 

u{x) = bia{x) sin + 62] , (5) 

where bi and 62 are complex constants. A straightforward substitution in Eq. @ leads 
to the two equations 

dla{x) + p^{x)a{x) = a{x) [dx4>{x)]'^ (6) 

where we can set = 1 without any loss of generality (since can be absorbed into 
a by redefining a a/a^/^), thus recovering Eq. (|). For obvious reasons, a and 4> 
are known respectively as the amplitude and phase functions. In terms of Ui and U2, 
it follows from standard results on Ermakov systems that the general solution for a is 
given by 



a{x) 



21 



+ 2I& 



21 

x) + 777:tU2(^) 777 ""1(^)^2(3^) 



1^2 



4/c 

w 



,1/2 



and the equation for the phase is readily integrated to give 

1 



arctan 



21 



2Ic' 



W^^ - 2Ic 

U2[X) 



+ arctan 2 Jc, 



(8) 



(9) 



where the integration constant is chosen so that = 0. Eq. (|) is an illustration 

of the nonlinear superposition principle 0. / and c are two constants, independent of 
X (J is the Ermakov, or Lewis- Ray- Reid, invariant). Note that the value of the phase 
function at S2, known as the accumulated phase, does not depend on the constants / 
and c for the eigenvalues Eq of the Sturm-Liouville problem, since 



{S2, E = Eo) 
(S2, E ^ Eo) 



Tin, 



arctan \2Ic[E)\ H vr; 



(10) 
(11) 



here n is an integer giving the number of nodes of U2- Note also that a{xi) is independent 
of c if xi is a zero of ui. 



2.2. Boundary conditions and energy dependence 

The boundary conditions for a and (f) are therefore incorporated through the parameters 
I{E) and c{E), which as indicated depend on the energy. Normalization of the 
eigenfunctions / of the Sturm-Liouville problem with vanishing boundary conditions 
at si and S2, yields 

r f{x)dx = IdE(t){s2, E = Eo). (12) 

Jsi 
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By choosing the eigenf unctions to be normahzed per unit energy increment, I becomes 
an energy-independent positive constant, and amphtude and phase functions depend on 
the single parameter c{E). 

We now introduce another solution g{x) of the Schrodinger Eq. (H), defined in 
terms of the solutions regular at si and S2 by 



g[x,c) 



21 



W 



cui{x) 



(13) 



1/2 



which fulfills W[mi, g] = 21 and gives a{x) = jj {ul{x) + g'^{x)) . This simply means 
that Ui and g lag 7t/2 out of phase, and that with the conventions of Sec. 2.1 (a = 1 
and 0(si) = 0), we have 



ui{x) = V 2/a(x, c) sin 0(x, c), (14) 
g{x,c) = V2Ia{x,c) cos (j){x,c). (15) 

We have emphasized the c-dependence of the different functions (though self-consistency 
requires it, it may be checked explicitly that ui does not depend on c). 



2.3. Oscillatory properties 

2.3.1. Oscillations of the amplitude Let ti{E) (t2{E)) be the inner (outer) turning 
point. For values beyond the turning points (when x < ti or x > t2) , we recast Eq. (P) 
as 

=/(x) -a"^(x), (16) 

where {(p;x) = dl(j)/dx(p — ^{d^cp / dx(py denotes the Schwartzian derivative. Since a 
is a positively defined quadratic form and p^(x) < beyond the turning points, the 
Schwartzian derivative of the phase is negative, and as it can be verified, if (0; a;) < on 
an interval then dx4> cannot have a positive local minimum on this interval, ie a cannot 
have a local maximum. By noting that a{x) +oo when x —* Si and x — > S2, we 
conclude that if a does not oscillate between the turning points, it will not oscillate on 
the whole interval ]si, S2[. Note that there is then a unique value xq, with ti < xq < ^2, 
such that d^a^Xo) = 0. 

However, as ui{x) and U2{x) oscillate for ti < x < t2, a will generally oscillate, by 
virtue of the nonlinear superposition principle Eq. (^, between the turning points, the 
local wavelength being half that of ui or U2. Nonetheless, at a specified energy, there 
may be infinitely many values of c giving a non-oscillating amplitude function||[ We 
give in the following paragraphs a sufficient (but not necessary) condition for c{E), that 
is the value of c and its explicit energy-dependence, because this value has remarkable 
properties related to classical quantities in the semiclassical limit, as will be seen in Sec. 
3. 

I For example, taking the second derivative of a, it is seen that any value of c enclosed between c±{x), 
with c±{x) — U2{x) [Wui{x)]~^ ± [2/1*1(2:)] ^ [2I/p{x) — uKx)]^^"^ and where x spans the interval 
between the turning points, will do. 
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2.3.2. Inverted phase accumulation Proceeding as in Sect. 2.1, we can define a phase 
function 0(a;) sucli tliat 0(s2) = 0, ie tlie pliase starts accumulating at $2 instead of Si. 
As can be easily seen, this amounts to exchange the roles of ui and M2; for instance, we 
now have U2 = 62«sin0, where (9^:0 = q;(x)~^ and 62 is a constant to be set below, a is 
given by 

with W['U2,'Ui] = —W and where we assumed 1 = 1 for simplicity. Generally a{x,c) 
and a{x, c) are very different functions (for example the c independent points of a are 
now located at the zeros of U2). Notwithstanding, it is apparent that a{x,c) = d{x,c) 
iff = = W-^ - (2/) -2 and c = -c. We shall set 



, , , , 



[17) 



Co{E) 



[W{E)r" - [2/] 



1-2 



1/2 



(18) 



To keep the quadratic form real, this implies that W"^ < 4P, condition to be assumed 
in the rest of the paper (this is not a problem in practice because a and are left 
unchanged by the transformations ui — > nui, W — >• kW, I —>■ k^/, c — >• c/k^, so the 
Wronskian can be conveniently rescaled). 

Let us now suppose that c = — c = ±Co. We then have a = a. 62 is found by 
evaluating W[mi, M2] at S2, which yields, by choosing a proper sign convention 62 = V^I- 
From Eqs. (|l)-(|ll), it follows that 

1/2 



sm I 



[TCo) 



W 



cos0(±Co)/2/ =F 



W~ 



{2iy 



sin 0(±Co) 



(19) 



thereby obtaining the relation between and 
only if a does not oscillate. 



which give the oscillations of ui and U2 



2.3.3. Auxiliary quadratic form on the unit circle a^{x) is a positive definite quadratic 
form. Labelling M the matrix of the coefficients, we have det M = and 
Tr M = 1/2/ + 2/c^ + 21 /W^. a'^{x) can be reduced to the canonical form 

a^(x, c) = Ai(c)f^(x, c) + \2{c)vl{x, c), (20) 

where Aj(c) are the eigenvalues of M (Ai > A2) and Vi{x,c) are the eigenvectors, 
normalized so that vf{x,c) + f |(x, c) = ul{x) + ul{x). We now introduce a quadratic 
form Q defined by 

Q{x, c) = Xi{c)wl{x, c) + X2{c)wl{x, c), (21) 

where wf{x,c) = vf{x,c) [^1(2;) +^2(2;)] As indicated, Q, as well as the Aj and Vi 
are c-dependent. Q oscillates between its maximum and minimum values, which are by 
construction given respectively by Ai and A2. 

We now set c = ±Co. Let xi {X2) label the points where ui (^2), vanishes; we then 
have Q{xi^2,^Co) = 21 /W^, so that between 2 zeros of Ui and U2, Q{x,±Co) has at 
least one extremum on the unit circle (for those points, the equality uKx) = ^2(2;) is 
fulfilled). Note that a{x, c = ±Co) goes through both the c-independent points of a{x, c) 
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and a{x, c), so if a oscillates, then there is an extremum of a between xi and X2, and the 
sign of d^a alternates between the consecutive zeros of ui and U2- This is illustrated on 
Fig. 1 for the specific case of the harmonic oscillator (to be discussed in details in Sec. 
5.4); the zeros of Ui and U2 are respectively shown as triangles and rectangles. Note also 
that the maxima of Q{co) correspond to the minima of Q{—Co), since Q{—Co) = Q{co)- 
Combining Eqs. ( PUD and (^Tj) and taking the derivative ^^^^(xi 2) as a function of d^Q, 
Q, Ui and U2, it can be seen indeed that for c = — Co the sign of d^a at 2 consecutive 
zeros, Xi and X2, alternates. However, for c = Co, the sign of dx<y between 2 consecutive 
arbitrary zeros of Ui and U2 does not change, and thus a(co) does not oscillate|. 

2.3.4- Oscillations of the accumulated phase We have explained in [§| why obtaining 
non-oscillatory functions is important when amplitude-phase methods are employed in 
scattering theory. The goal there is to extend energy-normalization, which for the 
eigenfunctions of the Sturm-Liouville problem is given by Eq. (|12D, to functions f{x) 
which converge at S2 but diverge at si (the phase-shifted or scattered wavefunctions). 
By combining the continuity equation for the probability density and L'Hopital's rule, 
improper energy normalization is defined by 

fix, E)dx = Wec + 2Ic\E)^ ' , (22) 

where r is a cut-off radius (and as above, I is assumed to be energy- independent). This 
normalization is of course arbitrary, since it is governed by c{E), but it conditions the 
energy-dependence of the different scattering parameters. In particular, the accumulated 
phase, which is unambiguously defined (Eq. ([To|) ) for the eigenfunctions of the Sturm- 
Liouville problem, crucially depends (Eq. (pT])) on the normalization when E is not an 
eigenvalue. 

More precisely, let us assume the eigenvalues of the Sturm-Liouville problem to be 
given by Eq = C,{n), where n is the number of zeros of the corresponding eigenfunction 
(thus, of U2), and ^{n) is an a-priori arbitrary, but monotonous function admitting a 
different iable inverse, n{E) = ^~^{E). The functional relation between Eq and integer 
values of n is thereby extended to any energy E lying between two eigenvalues, ie 
E = ^{n), n real. The energy-dependence for c{E) can now be chosen so as to extend 
the normalization of the eigenfunctions to non-integer values of n by equating Eq. (^) 
to IirdE^^E) (cf Eqs. and (0)), yielding 

c{E) = -jjCotnr\E). (23) 

Substituting in Eq. ( pT]) gives the following expression for the accumulated phase: 

0(S2, E^Eo) = nr\E) = 7rn{E). (24) 

§ We noted, however, that given the behaviour of the amphtude function at si and 82, there is 
necessarily a point between ti and t2 where dxCt vanishes. The argument sketched here rehes on 
the signs of the basis functions Ui and 1*2 and their derivatives when V{x) is monotonous on a full 
cycle of oscillation of the basis functions, and excludes the neighbourhood around the bottom of the 
potential, where dxCt changes sign (but dxCt does not vanish exactly at the bottom of the potential, as 
would be the case in the WKB approximation). 
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Thus the accumulated phase does not oscillate as a function of the energy (it is a simple 
straight line as a function of n) and the value of c given in Eq. (^3|) is the only value 
compatible with energy normalization leading to a non-oscillating accumulated phase 
function. We shall mention in Sec. 5.3 below the relation the specific form (^) has 
with the canonical action variable in classical mechanics. Note finally, that provided 
the basis functions Ui and U2 are redefined so that their Wronskian is proportional to 
21 simTn{E), Eq. (^) becomes a particular form of the more general Eq. (plSf ): with 
such a choice, the amplitude a is a non-oscillating function of x and the accumulated 
phase 0(52) is a non-oscillating function of E. 



3. Semi-classical Ermakov system 

3.1. Asymptotic solutions to the linear equation 

The approximate solutions to the one dimensional Schrodinger equation when h ^ are 
well-known from the asymptotic theory of ordinary linear differential equations It 
follows from Sec. 2.3.1 that it is sufficient to consider the solutions between the turning 
points (ie for real p{x)). Real solutions are of the form 



uix) = — , sm 



(25) 



± / p{x')dx' + 02 

/p{x) L 

where ai and 02 are constants. Tilded (~) quantities will henceforth denote asymptotic 
(semiclassical) functions when these are to be distinguished from the corresponding 
exact quantum solutions. It is well-known from Hamilton- Jacobi theory that 

S{x, E) = ± j p{x\ E)dx' + as, (26) 

where S{x, E) is known as the Hamilton- Jacobi characteristic function or reduced action: 
the characteristics in phase-space are made up of the points {x^d^S). 



3.2. Ermakov system 



By direct substitution of a general asymptotic solution into the Schrodinger equation, 
and by labeling a and the semiclassical amplitude and phase functions, we obtain a 
semiclassical Ermakov system 



+ 

a 



p^ix) + -— {S]x) 



u 







p^{x) 



+ 



{S;x) 



h^ {d^j 



(27) 
(28) 



where again the bracket ( ; x) denotes a Schwartzian derivative and we have, as for the 
usual amplitude-phase functions 5^ = a/da4> and thus d'^a/a = — {^0; /2. Eqs. (W^)- 
(p8|) are the semiclassical version of the quantum system given by Eqs. (![)-(§)• Eq. (^) 
is the modified Schrodinger equation exactly obeyed by the semiclassical wavefunctions, 
and Eq. (^) is the nonlinear equation fulfilled by the semiclassical amplitude function. 
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The passage from the exact (quantum) Ermakov system to the semiclassical one simply 
consists in a redefinition of the potential energy function, and is identical if {S; x) 
vanishes. Although generally {S;x) is nonzero (except for the free particle), {S;x) does 
tend to zero or to a finite value in the limit of high quantum numbers (eg the harmonic 
oscillator for the former, the centrifugal Coulomb potential for the latter). Only if this 
value is negligible compared to the other terms in the energy function does the high 
quantum numbers condition fit with the semiclassical limit. 



3.3. General solutions 

a and are given in terms of two independent functions Ui and U2 of Eq. (p^) 
by the same relations, Eqs. (H) and (H), as in the exact (quantum) case, with now 
tilded quantities. It is convenient, however, to set the tilded constants a, I and W 
equal to their quantum counterpart a, / and W. This is done by first noting that 
W[a sin 0, a cos 0] = a, which we then set equal to W[a sin 0, a cos 0], so a = a = 1. To 
preserve the Wronskians (cf Eqs. (P^-(P^), the semiclassical function ul, of the form 
given by Eq. (^^), is thus set as 



= y^sin5(x), (29) 

where we have implicitly included 02 in the reduced action so that ul is the asymptotic 
approximation to Ui in the neighbourhood of an arbitrary x lying between the turning 
points. An independent solution with Wronskian W[ui,U2] = W is then obtained 
under the form cos(S'(a;) + b)/Jp{x) as 



U2[X) = J — — KCOS 
V P{^) 



w 

S{x) + arccos 



2M ■ 

where we have introduced the scaling factor k to keep all quantities real. In what 
follows, we shall set k = 1, which is tantamount to rescaling U2 (a similar rescaling was 
performed in the quantum case, see below Eq. ([T8|)). 



3.4. Non-oscillating solutions 

Substituting Eqs. (p9D -(^) in the expression for a readily yields 

'aI^W-^ sin^ S{x) + (cos5(x) - 2/csin^(x)) 

1/2 



a 



m-i 
1 

p{x) 



X 



COS S{x) - 2/ c + W'^ 4 - 



sin S{x] 



(31) 



which is a highly oscillatory function for an arbitrary value of the parameter c. However, 

1 /2 

it may be noted by inspection that for c = — [W~'^ — (2/)^^] , the oscillating terms 
are cancelled out. Remark that this is the same expression that was labelled Co in the 
quantum case (Eq. (|18D). Reestabhshing h, the amplitude now reads a^(co) = h/p{x), 
which given our assumptions is a non-oscillating function of x. Identical substitutions 
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may be done for 0, from which it follows that the semiclassical phase function is highly 
oscillatory for an arbitrary value of c except if c = Cq, and in that case, 

^{x,Co) = Six)/n. (32) 

In short, non-oscillating functions are obtained for a unique value of the parameter 



c, for which the semiclassical quantities match their classical counterpart {\Jh/p{x) 
and S{x)/h are respectively the classical probability amplitude and phase functions). 
Writing Eq. (pSf) as 

y = /(a:) - (33) 

this means that each side of the nonlinear equation of the semiclassical Ermakov system 
vanishes independently. Note also that p'^{x) + h"^ {S; x) /2 = (tan(p{c);x^ /2 (this is 
established by using the Mobius invariance of the Schwartzian derivative and establishing 
a c-dependent linear transformation relating tan0(x, c) to tan5'(x)). 



4. h expansions 

The link between the solutions of the quantum and semiclassical Ermakov systems is 
done by employing a formal h expansion of the amplitude function. As in the previous 
section, we assume / and W to be identical in both the quantum and the semiclassical 
case. It is then straightforward to show that to each order in h (as well as to infinite 
order), non-oscillatory functions are obtained. 



4-1- Series expansions 

The formal asymptotic solution to the Schrodinger equation for small values of the 
parameter h is usually done by transforming it to the Riccati form and then obtaining a 
recurrence relation between complex function of order j and the functions of lower order 
p!0| . Here we proceed slightly differently, because we want the relations between the 
amplitude and the phase to be verified to each order. We look for a generic solution of 
Eq. (^ under the form u{x) = a{x) exp i f (x) / h where a{x) and f{x) are real functions 
admitting the series expansions 



j=0 3=0 

Substitution into the Schrodinger equation gives aj{x) = 
following recurrence relations for even j, j > 2: 



0, fj{x) 



(34) 

for odd j and the 



dxfj{x) 



aj{x) 



dla 



j-2 j-m 

j—2 ^ ^ '^m ^ ^ 9xfn9xfj—m—n 
m=2 n=0 



i-2 

E 

n=0 



/ (2dyandyfj_n + cindyfj-n) / \j dyf^dy 



(35) 
(36) 
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with /o(a;) = ±S{x), ao(x) = ±dxfo , and where the constants bj appearing in the 
solution to the homogeneous equations for aj{x) are all set to for j > 2. The Wronskian 
relations are then preserved to each order in h , that is 

W[a(x) sin /(x),a(a;) cos /(x)] = al{x)dr,fo{x) = dn,S{x) /p{x) = 1. (37) 

Between the turning points, the ± branches are combined to yield real oscillatory 
functions. The formal expansion for Ui (assuming again the adequate integration 
constant to be included in /o) is then given by 



uAx) 



oo / oo 

YlY^aj{x)W sm[Y^U{x)W 

j=0 



(3^ 



\i=0 



The expansion for the function lagging 7r/2 out of phase is trivially obtained by using 
the cos function; from Eqs. (|T3D to (p!5| ) the formal expansion for U2 is then found as 



U2[X 



21 



W cos (Y,Mx)h' 

\i=0 



4 - WyP 



1/2 



sin iJ2fi(^W 

\i=0 



(39) 



To first order these functions coincide by construction with the semiclassical 
wavefunctions ui and U2- 



4-2. Amplitude and phase expansions 



The h expansions for the amplitude and phase functions are obtained by combining the 
nonlinear superposition principle [Eqs. (|^)-(^)] with the formal series expansions for ui 
and U2- The expansions may be done to finite or infinite order. In the first case, the 
functions are Taylor expanded around ^ = after separating the classical terms oq and 
fo- The infinite order case is analogous to the first order case treated in Sec. 3.4. For 
example substituting the series expansions in Eq. (|) gives the expression (mod vr) of 
the phase which can be simplified as 



ix,c) = arccot < cot 



'S{x) 

h 



i-1 



2Ic{E) + IW-^ A-W'^/r 



1/2 



.(40) 



The amplitude function may be obtained by deriving this last equation, keeping in mind 
the Wronskian relations (^), yielding an expression involving the sines and cosines of 
the expression between parentheses in Eq. (^Of ). In both cases, the highly oscillatory 
terms are cancelled by setting c = Co (in Eq. (|^) for example, the term between square 
brackets then vanishes). 



5. Discussion 



5.1. General Remarks 

We have thus seen that same value of c gives rise to non-oscillating functions both in 
the semiclassical and quantum cases. This is not surprising, if the similarities between 
Eqs. (1^), (^) and ( |19|) on the one hand, together with the nonlinear superposition 
principle in both the quantum and semiclassical Ermakov systems on the other hand 
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are considered: the same functional relation gives the amplitude-phase functions in 
both cases. We note however that the simple arguments we have given above, although 
physically appealing because of the direct connection to standard classical quantities, are 
liable to a more rigorous treatment. For example the series expansion obtained for the 
amplitude function does not necessarily converge, and even its asymptotic properties 
as a function of the parameter c deserve a more thorough investigation. Actually, 
Lewis had made the same remarks when studying the adiabatic invariant series of the 



time-dependent classical harmonic oscillator in powers of an adiabatic parameter e [11 
which is defined by the same uncoupled Ermakov system as the one studied in this 
paper. Lewis's work was the first (albeit implicit) application to Ermakov systems 
of Kruskal's asymptotic theory of Hamiltonian systems extending the study of the 



adiabatic invariants beyond the first order in e ||T^ . The transposition of these theories 
to the present problem is not straightforward because we lack here the Hamiltonian 
formalism on which these theories are based (eg, the integral invariants that appear 
in Kruskal's theory would have here a rather obscure interpretation). Nonetheless, the 
present results on the oscillatory properties can be directly transposed to any uncoupled 
Ermakov system depending on an 'adiabatic' parameter e (which in the present context 
corresponds of course to the Planck constant). 

5.2. Scattering basis functions 



Previously to the work of Fano et al [|l^, the use of amplitude-phase methods in 



scattering theory was limited to the high kinetic-energy limit, and the derivatives of 
the phase above first order neglected, thereby effectively restricting the treatment 
from the start to the standard WKB approximation (see eg, Ch. 4.3 of [1^). The 



more recent application of these methods to define a pair of basis functions for phase- 
shifted wavefunctions in a potential (as outlined above), relies on numerical treatments 
to minimize the oscillations; these treatments are preferred even in the cases for which 
an analytic pair of basis functions is known, such as the Whittaker functions for 
the centrifugal Coulomb problem or the parabolic cylinder functions for the harmonic 
oscillator. It is interesting to note that the approach suggested in Sec. 2.3 yields, for 
an arbitrary potential with a single minimum, the same relations that are known to be 
valid (and non-oscillating) for the analytic functions in the mentioned special cases (eg 
in the Coulomb case where the effective quantum number u is defined by z/ = ^"^(i?) +/, 
the accumulated phase obtained with the Whittaker functions is (t){oo) = 7r(z/ — /), and 
improper normalization follows the normalization of the eigenfunctions by normalizing 
to i^^/2, independently of being real or an integer [|15|] ). 

^From a formal standpoint, defining a specific basis of functions is equivalent to 
defining the Green's function of the scattering process in the asymptotic field. This is 
the Green's function that appears in the Lippmann-Schwinger equations and through 
which the collision operators are defined. This is why the collision operators depend on 
the parameter c. Though at first sight this may appear as an unexpected feature, it must 
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be remembered that the exphcit inclusion of closed channels leads to a modification of 



the usual Green's functions through a term depending on the accumulated phase |16 



5.3. Semiclassical limit and classical quantities 

Standard semiclassical physics is usually not concerned by the quantum to classical 
limiting procedure, because the strategy there is to start from classical quantities at 
the outset and then proceed to quantization. However, we have seen in Sec. 3 that 
the classical reduced action is only one of the many phase functions that are obtained 
in the semiclassical limit (namely the non-oscillating one). This type of problem is 
frequent in 'classicalization' procedures: a particular, often arbitrary choice has to be 
made to recover classical quantities (action, Liouville equation, etc.) A recent example 
is given in Ref. [0, where the passage from Hilbert-space to classical- phase space 



operators involves particular choices for the parameters in order to recover the classical 
dynamics. In other works, this requirement takes the form of an additional ad-hoc 
condition usually termed as the 'correspondence principle' (eg in Ref. where 
a quantum version of the Hamilton- Jacobi equation is given an additional boundary 
condition Pguantum{x, E) Pdassicai{x,E) whcu ^ ^ and E is fixed)|]]. Though much 
less general, our treatment is more transparent in that the continuation of the classical 
reduced action in the quantum domain is readily identified: it is the non-oscillating 
phase function (and the continuation of the classical probability amplitude is the non- 
oscillating amplitude function). 

Other aspects of the classical-quantum correspondence for a classically integrable 
and separable system deserve to be mentioned. The accumulated phase (Eqs. ([TO|)-p!lD) 
is seen to be directly related to the line integral around a closed loop of : 

f[dMx,c)]dx = 2(l){s2,c). (41) 

Contrary to EBK (torus) quantization, there are no caustics when dealing with exact 
quantization, and the quantization condition reads / [dx(f){x, c)] dx = 27m, ie, n is an 
integer and the line integral does not depend on the particular value of c. Here we 
understand by "exact quantization" the quantization of the exact quantum phase, and 
not the exact WKB quantization of non-solvable potentials, as employed by Voros [l^ . 



Note however that when c = Co, the (unquantized) integral (|41|) reads, according to Eq. 



[5x0(x, Co)] dx = 27m{E). (42) 

Not only does c = Cq preserve for non-integer real numbers the functional relation valid 
for exact quantization, but it gives a parameterization of the quantum equivalent of the 
canonical action variable. 

II Other authors crudely suppress the ?i— dependent terms in selected equations where this suppression 
leads to classical relations (for example this would be done in Eq. (Q) by giving the amplitude squared 
the dimensions of a classical quantity, getting thereby rid of d^a, but would not be done in Eq. (^), 
which does not support an obvious classical interpretation). This procedure has often been criticized 
because the functions and the parameters appearing in the equations depend on h. 



Non-oscillating solutions to uncoupled Ermakov systems and the semiclassical limit 14 



In a similar vein, the period of motion T is given in Hamilton- Jacobi theory 
by taking the energy derivative of the reduced action along the closed loop. The 
transposition to the phase function in the quantum case would imply taking the 
energy derivative of Eq. (0), which by Eqs. ( ]Tl| ) and (|22| ) is proportional to the 



normalization. Again, the normalization depends on c{E) and to take dE(f){c) as the 
time parameterization doesn't appear to make much sense unless c = Cq, since any other 
value would lead to an oscillating function, which would further not collapse to the 
classical period in the h ^ limit. 

5.4- Example 

We illustrate the properties mentioned above on the harmonic oscillator, a paradigm 
both in the Hamilton- Jacobi formulation of classical mechanics and in the 

semiclassical theory of bound states The reduced action S{x,E) of Eq. (^) is 

readily obtained, from which it follows that 

d^S = (2mE - mWx^) (43) 

J = ^ [d^S] dx = 27tE/uj, (44) 

where m and u are the mass and frequency of the oscillator; we slightly depart from 
usual conventions and define J to be the canonical action variable. The period is 
recovered as T = OeJ- The standard WKB solutions between the caustics, Eqs. 
(p9|)-(|30|), are obtained from these classical quantities. Semiclassical quantization must 
take into account the singularities at the turning points, from which it follows that 
Eq = hu (n + 1/2) , with n being an integer. 

Quantum mechanically, the eigenvalues are given by ^(n) = hu {n + 1/2) when n is 
an integer that counts the zeros of the eigenfunction, but this relation can be extended for 
any real value of n, as discussed in Sec. 2.3.4. We then write as n{E) = E/huj — 1/2. 
The derivative of the exact (quantum) phase is given by Eqs. (0)-(|§) as 



„ , , . muj r / 1 2\ 2/ \ 2 J 2/ \ 4/c , , / n"' ^ 

d^(f){x, c) = — ( — + 2/c j u^{x) + ^u^ix) - —ui{x)u2{x) 



Fig. 2 compares the classical conjugate momentum as given by Eq. (|43| ) with two exact 
phases hdx4>{x, c) obtained from Eq. ( ^3]) by using numerical solutions of the Schrodinger 
equation ui and U2 with respective vanishing boundary conditions at x = — oo and 
X = +00. One of the curves is for an arbitrary value of c, the other corresponds to 
c = Co- It may be seen that even for a moderate excitation (n ^ 12), the non-oscillating 
solution with c = Cq can barely be distinguished from the classical momentum at the 
same energy, except near the turning points (the reason is that {S; x) is negligible, hence 
the exact quantum phase, solution of Eq. (^ tends to the semiclassical quantum phase 
of Eq. (|28|) , which is simply S{x) when c = Co). 

The analog of the canonical action J appears as the line integral (^T]), which 
generally depends on c except when n is an integer in which case / [hdx4>{x, c)] dx = 
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2T\hn. Thus for the quantized energies, the quantal hne integral differs from J by the 
Maslov index. However, if c = Cq, we have for any E [Eq. (^21)] 



j> [hd^(t){x, Co)] dx = 2'kE/uj - (46) 

which is the classical result with an action correction coming from the Maslov index. 
Note that taking the energy derivative of Eq. (l4l|) crucially depends, for any energy E 
(including the eigenvalues) on the energy-dependence c{E). Only the energy dependence 
given by the relation (p3|) above, which yields Eq. (|46|) in this case, renders the usual 
relation for the period, and more generally follows the classical time parameterization for 
conservative systems in Hamilton- Jacobi theory. Thus, the parameter c, which appears 
free within quantum mechanics and has no classical counterpart, must be constrained if 
the usual classical relations for the oscillator are to be extended to quantum amplitude 
and phase functions. 

The usual interpretation of the invariant within Ermakov systems hinges on the 
use of an original Hamiltonian or Lagrangian, from which the Ermakov equations are 
derived. The invariant is then associated, by means of Noether's theorem, with the 
conserved quantity of an auxiliary motion [|, ^ . Such an interpretation is of course not 
available here, where the Ermakov equations are used in a quantum-mechanical context. 
The invariant / is employed in this context only to define the normalization, through 
Eq. (12). However, for the choice further interpretative step may be taken. 



since then the term on the right hand-side in Eqs. ([T^ ) or (|2^) gives 
h f 

—IdEj[hd:^4>{x,Co)]dx. (47) 
For the harmonic oscillator, we have by Eq. (^) and by adopting unity normalization 



{At 



hn 

For other systems, / usually depends on n{E) except if the wavefunctions are energy- 
normalized. Elementary manipulations yield the more general form 

-1 



2n 



X{x, E)dx 



(49) 



where A(x, E) is the local de Broglie wavelength and the range of integration is restricted 
to the classical domain between the turning points. 



6. Conclusion 



Previous interest in amplitude-phase methods led us to investigate in this work the 
oscillatory properties of the nonlinear equation of uncoupled Ermakov systems. It was 
shown that non-oscillating amplitude-phase functions in the space and energy variables 
have a particular feature in the semiclassical limit: they yield classical quantities. 
We have seen that although standard quantum-mechanical quantities, such as the 
wavefunctions or the eigenvalues, are insensitive to the value of the parameter c and its 
energy dependence, there is a unique value of c which appears as connecting quantum 
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amplitude and phase functions to their classical counterpart. Only in this case can 
'quantum characteristics' {x, dxcf)) having a sense, and collapsing to (x, d^S) when Ti ^ Q 
be defined. We insist again that from the point of view of quantum mechanics, even in 
the semiclassical limit, this need not be the case: any amplitude and phase functions 
obeying Eq. (^) will yield correct semiclassical wavef unctions. A full study of these 
aspects on specific physical systems will be given elsewhere. Further links with current 
work on Ermakov systems may lead to a better appreciation, as well as to an extension, 
of the formalism. 
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Figure captions 



Figure 1. Q{co) (broken line, left scale), a(co) (non-oscillating black line) and a(— Co) 
(oscillating grey line) [right scale] are plotted for a harmonic oscillator (atomic units 
and oj = 1, n{E) = 4.4, see Sec. 5.4 for the definitions). All amplitude functions a{c) 
for any c go through the points xi (triangles), and any function a{c) goes through the 
points X2 (rectangles). Only the 2 functions a(±Co) go through both the points xi and 
X2- Since Q oscillates and Q{x = a;i, 2, ±Co) = 21 /W'^ is constant, dxQ has opposite 
signs at Xi and X2, as may be seen on the figure. This is also the case for a(— Co) (grey 
line) which therefore oscillates, but not for a(co). 



Figure 2. The positive branch of the classical canonical momentum p{x) for a 
harmonic oscillator (w = 1, atomic units) is plotted (grey line) vs the "quantal 
momentum" dx^{x,c) i) for an arbitray value of c (dashed line) and ii) for the non- 
oscillating value c = Co (black line). 



